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The hydrodynamics of the superfluid interior of a neutron star--and 

associated relaxation phenomena presumably susceptible to observation~are 

quite sensitive to the isotropic energy gap A. of the neutron matter of 

k F 

the crustal layers [1]. The most quantitative evaluation of A. is the 

r 

variational calculation of Yang and Clark [2] based on a wave function 
incorporating short-range Jastrow correlations as well as longer-range 
pairing or BCS correlations. For pure neutron matter and an adequately 

realistic two-neutron potential, these authors find that A,, peaks at 

-1 F 
k c = 0.72 fm , with a value A. = 2.45-MeV. 
r k f 

It has been argued by Pines [3] that this kind of evaluation is likely 
to yield a result for the energy gap which is on the low side, because the 
polarizability of the neutron medium has been essentially neglected. Ac- 
cording to an estimate by Pines and Pethick, polarization of the medium 
tends to enhance a "bare" attractive interaction between two neutrons 
with wave vectors k, -k near the Fermi surface. The enhancement factor 
was found to be (1 + F where F Q is the leading Landau Fermi-liquid 
parameter appearing in the Legendre expansion of the spin-symmetric part 
of the interaction between two quasiparticles on the Fermi surface. With 
F q ^ -0.7, a substantial amplification of the pairing matrix elements would 
result, and the energy gap aid condensation energy, which are extremely 
sensitive functions of these matrix elements, would increase dramatically. 
Indeed, these considerations raise the intriguing possibility that polariza 
tion and pairing may conspire to bring about a first-order phase transition 
in low-density neutron matter or even a bound state, metastable or stable. 

In this note we shall take a closer look at the singlet-state quasi- 
particle interaction in neutron matter. Our considerations will be based 
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on the results of a detailed evaluation of the Landau Fermi -liquid para- 
meters for pure neutron matter [4] including polarization effects, i.e., 
including the interaction induced by exchange of density and spin-density 
excitations. It will be concluded that, owing to the spin-dependence of 
the quasiparticle interaction, and ultimately the balance of attraction, 
repulsion and spin-dependence in the fundamental two-neutron interaction, 
polarization actually works to suppress rather than to enhance the pairing 
matrix elements. 

In the theory of Babu and Brown [5] the Landau quasiparticle inter- 
action energy is approximated as 

^ (&1 »! S 2 * *~2 ^ ~ ^d^~l ’~2*~1 *~2^ M~1*~2’ 2 1 ’~ 2 ^ * 

where the "direct" part ^ is obtained by functional differentiation of 
the lowest-order Brueckner approximation to the ground-state energy and 
the "induced" part of ^ aiises from exchange of density and spin-density 
excitations. The induced interaction between quasiparticles with wave 
vectors k^,]^, ln individual spin states s^.Sg, is given by 


S 1 S 2, 


S-jOt 


/ fin 


V ‘a r k 2 ) = - L < ' %,£) 
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Bs, 


fiu 


p'.B 


i (p‘,k 2 ). 


( 2 ) 


q»0 


and represented graphically in fig. 1, where the blob stands for the 
response function (fin /fiu , J n . 3abu*and Brown furnish a transport 

P*& p >P QiU 

equation for (fin /fiu , „) in terms of the full interaction £ (p»P')* 

B ,a B ~ ~ 

(N.B. fin p a is the change in the quasiparticle occupation number n p q due 

to a weak external potential fiu . Q (q,a)), where q,u> is the four- vector 

p p - 

of the induced disturbance.) Strictly, (2) applies to the long-wavelength 
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limit q = JjCj-kgl * 0* an extrapolation from this limit is performed as 
described in refs. [4,6], permitting (2) (with use of (1) and the trans- 
port equation) to be solvedfor the induced interaction ^ . Certain 
diagrams must be omitted from the direct part of & so as to avoid double 
counting. 

The spin dependence of the quasiparticle interaction takes the form 

£ (ij'j »!$2 *2^ (IS-j »iS2 ^ ^ 9 ( !S^ ® l$2 ^ ~ "l *~2 * 

and similarly for and The spin- symmetric part f and the spin- 
antisymmetric part g may be expanded in Legendre polynomials of 
x e cos(k-j,k 2 ): 

f(k r k 2 ) = £ f t p iM * = £ g i?i (*) • (4) 

With ki = k 2 =kp, the Landau parameters f^, g^ depend only on the density. 

Keeping just the first two terms in these Legendre expansions, Backman 
£t §1. [4] derived a useful approximate expression for i(. . Numerically it 
turns out (in their calculation for the Reid potential) that the parameters 
fj and g.j are small compared to f Q and g Q in the density range of most 
interest for isotropic superfluidity (kp = 0.4-0.8 fm"^). It is therefore 
a good approximation to drop terms in and g^, and work with this 
simplified version of the expression of ref. [4]: 
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Here we use an arrow notation for the individual spin states s^, 
s^ . The capital F*s and G’s are dimensionless , being obtained 
from the lower case f’s and g's by multiplication with the den- 

0 9 

sity-of-states factor N(0)=m*k^/7itr where m* is the quasiparticle 
effective mass. Finally, U(q) is the Lindhard function 



The contribution of the induced interaction to the (dimension 
less)singlet quasiparticle interaction is given by the spin- 
symmetric part of N(0)^ minus 3 times its spin-antisymmetric 
part (since 6.-* -3) , i. e., by 

_ I 

2 

- 3 

(v/e note that" at these densities D and higher even waves are 

unimportant./ Taking numerical values for F n and G n from ref. 4 

-1 1 

at kj l ry -'0.5 fm , we find, to a very good approximation, F^( s Q ) 

= 1. 2 TJ(q). This is a smooth function of q on the individual 
£ 0, 2k F J . Performing an angle average (a qdq/2k ? ^ integration 
over the stated interval), we arrive at ( s Q ) = 1.2 IT ( 1 .5kp)=0.9. 
This result is almost independent of density in the aforementioned 
k 7 range. 

If the second term in square brackets in (7) were negligible 

* 

compared to the first, cur result would conform with the original 

estimate - |V? |/(l+F ) of Pines and Pethick, if U(q) were set 
o o 

unity(not far off) and the^bare" interaction V identified with F q . 

However, with G ^ 0.8 and F~-0.3 according to (4) , the second 
’ o o 
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term clearly dominates, reversing the sign of the polarization- 
induced interaction relative to the Pines-Pethick estimate. 

vVe are indebted to C, Pethick for the following clarifying 
remarks: The main reason the present conclusion differs from that 
of ref. [7] is that the calculated equation of state of neutron 
matter has become less unstable against density fluctuations 
v'F has grown larger) since the time of the Pines-Pethick estimate 
The latter was made including only the density fluctuation channel, 
since it would be dominant over the spin-1 exchange channel for 
F q — >-1. ( This is very similar to the paramagnon model for 
liquid^He, except that there & — >-1.) Of course, if F q were close 
to -1, polarization effects would still enhance the singlet quasi- 
particle interaction even if the spin-dependence of the quasi- 
particle interaction were taken into account. 


For the direct part of the quasirariicle interaction we may 
take the renormalized form [6} 


S 1 S 2, 


i (JSi*JS 2 > = 0-2fC><k 1 s i s ? itikiSi,k 0 s 0 -ex> , (8) 

d t 

where t is the Brueckner reaction operator defined in terms of 


the usual choice of single-particle spectrum having a jump at the 
Fermi surface [ 7 ) and K is the wound parameter. At k^~* 0 . 5 fm 1 , 
the direct part of the(dimensionless) singlet quasiparticle inter- 
action is accordingly estimated to be F^( s q ) = -3.5. The net 
effect of the polarization contribution is thus to suppress the 
singlet interaction by a factor^0.7-0.8 relative to f.,(^e 0 ). 

The quantity F d ( 1 s Q ) depends more strongly on density than F i ( 1 s o ); 
hence the suppression factor {3 will show some appreciable density 


dependence . 
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</hat effect will ’the associated supression of the pairing 
matrix elements have on the energy gap and condensation energy 
of neutron matter? A simple-minded answer may be based on the 
weak-coupling formula (8) = (4fi 2 kp/m*) exp(-1 /N(0) V) . This 

formula is used once, inserting the 4 kp result of Yang Clark 

(?.') and their m* value, to determine a "bare” V, then again, with 
V — |3V, to calculate the suprcssed gap. At kp=0.6 fm”\ we find, 
taking ^=0.74, that the energy gap is cut down from 2.24 Mev to 
0.69 MeV. Correspondingly, the weak-coupling approximation predicts 
that the condensation energy E is suppressed! from 0.301 MeV) by 

C 

an order of magnitude. An elaborate evaluation of gap and conden- 
sation energy using the full method of Yang and Clark, with the 


pairing matrix elements P^ of that approach replaced by £P^ 


(but no other modifications), yields the suppressed values^. - 

k F 

0.64 MeV and E = 0.03 Mev. This evaluation rests on the unjustified 
c ° 


but at first sight not implausible identification of the effective 


interaction <1 2 Iw^l 1 2— 2 of the Yang-Clark method with the 

direct quasiparticle interaction. It is important to remember , however , 


that the calculation of ref. C2J is variational in nature, 


dealing with the expectation value of the raw neutron-matter 
Hamiltonian. (The three-body and higher-order cluster contributions 


to the Hamiltonian expectation value, not treated, are almost 
certainly negligible at these low densities.) Therefore the gain 


of energy E which Yang and. Clark obtain with their Jastrov-BCS 

V 

trial wave function, over the simple Justrow-Fermi gas energy 
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evaluation, is surely genuine. However, it could be that the 
Jastrow wave function used for the normal ground state does not 
adequately incorporate the effects of low-lying virtual excita- 
tions( especially those corresponding to polarization of the 
medium). It could be that an improved superstate trial wave 
function, incorporating short-range correlations, polarization 
effects, and pairing correlations, would lead to essentially the 
same or a somewhat increased energy gain, but with a substantial 
reduction in the optimal energy gap A, . ’.Ye are currently looking 

Kj, 

into this possibility. 

Evidently a workable first-principles theory of pairing in 
the presence of both short-range correlations and polarization 
is needed. 


# 
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Figure Caption 


Fig. 1. Induced interaction (spin labels omitted for simplicity). 
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